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THE HEUN EQUATION AND THE 
CALOGERO-MOSER-SUTHERLAND SYSTEM II: 
PERTURBATION AND ALGEBRAIC SOLUTION 


KOUICHI TAKEMURA 


Abstract. We apply a method of perturbation for the BC\ Inozemtsev model 
from the trigonometric model and show the holomorphy of perturbation. Con¬ 
sequently, the convergence of eigenvalues and eigenfuncions which are ex¬ 
pressed as formal power series is proved. We investigate also the relationship 
between L 2 space and some finite dimensional space of elliptic functions. 


1. Introduction 


In this paper, we report some properties for eigenvalues and eigenfuncions of the 
BC\ Inozemtsev model. Consequently, we obtain results on Heun function. 

The BC\ Inozemtsev model is a one-particle model of quantum mechanics whose 
Hamiltonian is 



(i.i) 


where p(x) is the Weierstrass p-function with periods (l,r), wo = 0, uj\ = 1/2, 
W 2 = (1 + t)/ 2, lo 3 = t/2 are half-periods, and l t , (i = 0,1,2,3) are coupling 
constants. This model is sometimes called the BC\ elliptic Inozemtsev model, 
because the potential is described by use of elliptic functions. 

There are two evidences which ensure the importance of the BC\ Inozemtsev 
model. The first one is equivalence to Heun equation, which will be explained 
in section [3 The other one is that BC\ quantum Inozemtsev model is a special 
(TV = 1 ) case of BCn Inozemtsev model (3, which is a generic integrable quantum 
system with Bn symmetry. In fact, classification of integrable quantum systems 
with Bn symmetry was done by Ochiai, Oshima, and Sekiguchi [5j, and it was 
shown that integrable quantum system with Bn symmetry is BCn Inozemtsev 
model or its degenerate one. It is known that BCn Inozemtsev system contains 
the well-known Calogero-Moser-Sutherland system with Bn symmetry as a special 
case. 

In this paper, we try to obtain physical eigenfunctions and eigenvalues of the BC\ 
Inozemtsev model, and investigate their properties. Here a “physical” eigenfunction 
means that it is contained in an appropriate Hilbert space, which is often a space 
of square-integrable ( L 2 ) functions. Note that roughly speaking the “physical” 
eigenfunction corresponds to the Heun function of Heun equation. 
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Applying a method of perturbation is a possible approach to this problem, which 
was done in El El for the Calogero-Moser-Sutherland system of type Am- Now we 
explain this method shortly. Elliptic functions have a period r. By a trigonometric 
limit p = exp(7r-\/—lr) —> 0, the Hamiltonian of the BC\ elliptic Inozemtsev model 
tends to the Hamiltonian of the BC\ Calogero-Moser-Sutherland model, and it 
is known that eigenvalues and eigenstates of the BC\ Calogero-Moser-Sutherland 
model are obtained explicitly by use of Jacobi polynomials. 

Based on eigenstates for the case p = 0, we can obtain eigenvalues and eigen¬ 
states for the BC\ elliptic Inozemtsev model (p ^ 0) as formal power series in p. 
This procedure is sometimes called an algorithm of perturbation (see section HOI . 
Generally speaking, convergence of the formal power series obtained by perturba¬ 
tion is not guaranteed a priori, but for the case of BC\ elliptic Inozemtsev model, 
the convergence radius of the formal power series in p is shown to be non-zero (see 
Corollary Id. 81) . and this perturbation is holomorphic. As a result, real-holomorphy 
of eigenvalues in p and completeness of eigenfunctions is proved. Note that a partial 
result was obtained in part im by applying Bethe Ansatz. 

There is another method to investigate eigenvalues and eigenstates of the BCi 
Inozemtsev model. If the coupling constants Iq , h. h, h satisfy some equation, 
the Hamiltonian H (see 111. H i preserves a finite dimensional space of doubly pe¬ 
riodic functions which is related to the quasi-exact solvability El El On a finite 
dimensional space, eigenvalues are calculated by solving the characteristic equation, 
which is an algebraic equation, and eigenfunctions are obtained by solving linear 
equations. In this sense, eigenvalues on a finite dimensional space are “algebraic”, 
and eigenvalues and eigenfunctions on a finite dimensional space would be more 
explicit than ones on an infinite dimensional Hilbert space. 

In this paper, we also investigate relationship between Hilbert spaces ( L 2 spaces) 
and invariant spaces of doubly periodic functions with respect to the action of the 
Hamiltonian H. In some cases, a finite dimensional invariant space becomes a 
subspace of the Hilbert space. Then it is shown under some assumption that the 
set of eigenvalues on the finite dimensional invariant space coincides with the set 
of small eigenvalues from the bottom on the Hilbert space. 

This paper is organized as follows. In section 0 the relationship between the 
Heun equation and the BC\ Inozemtsev system is clarified. Next we consider a 
trigonometric limit and review that eigenstates for the trigonometric model are 
given by hypergeometric (Jacobi) polynomials. We also explain how to apply an 
algorithm of perturbation in order to obtain formal eigenvalues and eigenfunctions 
for the BC\ elliptic Inozemtsev model. In section 0 we prove convergence of the 
algorithm of perturbation by applying Kato-Rellich theory. We also obtain several 
results related to Kato-Rellich theory. Although holomorphy of the eigenfunctions 
v m {x,p) in p as elements of L 2 space is shown by Kato-Rellich theory, conver¬ 
gence of the eigenfunctions for each x is not assured immediately. In section 0 we 
show uniform convergence and holomorphy of the eigenfunctions v m {x,p) for x on 
compact sets. 

In sectional finite dimensional invariant subspaces of doubly periodic functions 
are investigated and relationship to the Hilbert space (A 2 space) is discussed. In 
section 15.31 we focus on the case lo,h,h,h € Z>o- In section El examples are 
presented to illustrate results of this paper (especially section 0. In section 0 we 
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give some comments. In section 0 some propositions are proved and definitions 
and properties of elliptic functions are provided. 

We note that some results of this paper are generalized to the case of the BCn 
Inozemtsev model (see m- 


2. Heun equation, trigonometric limit and algorithm of perturbation 

2.1. Heun equation and Inozemtsev system. It is known that the Heun equa¬ 
tion admits an expression in terms of elliptic functions and this expression is closely 
related to the BC\ Inozemtsev system da ©Emm. In this subsection, we will ex¬ 
plain this. 

Let us recall the Hamiltonian of the BC\ Inozemtsev model 

d 2 3 

H = — ——~b FI + l)p(<£ + Wj), (2-1) 

2=0 

where p(x) is the Weierstrass p-function with periods (1 ,t), wo = 0, uj\ = 1/2, 
W 2 = (1 + t)/2, lo 3 = r/2 are half-periods, and U (i = 0,1,2,3) are coupling 
constants. Assume that the imaginary part of r is positive. Set 


ei = p(ui) (i = 

lp + 1 

&(w) = w 2 (w 


1,2,3), 


a = 


e 2 ~ e 3 
ei - e 3 : 


Ii+l io + l 


Note that a is nothing but the elliptic modular function A (r). We change a variable 

by 

ei — e 3 


p{x) - e 3 ' 


( 2 . 2 ) 


Then 


i ~ ( d 2 

o Ho = —4(ei — e 3 )<! w(w — 1 )(aw — 1)(-— ) 

l K dw' 


+ 


1 ^ 2/q “h 3 2 1\ -\- 3 ^ 01 ( 2/2 “h 3) ^ d ^ 


): 


aafiw + g j, 


(2.3) 


_ lp+ll+^2+^3+4 


where q — (pp JA_ 0 Zj(Zj + 1 ) — a(Zo + h + 2) 2 — (Iq + Zi + 2 ) 2 ), a — 2 

0 — ^ 0 +^l+^ 2 — ^3+3 

Let f(x) be an eigenfunction of H with an eigenvalue E, i.e., 

j 2 3 

(H - E)f(x) = j ~2 + '^2 l i{h + 1 )p{x + u>i) - e'J f(x) = 0. 


(2.4) 


2=0 


From (EHl) and (IQ) . we obtain 


(, d ,2 /^o + § h + | h + § x d 


—)-h 

1 / 27,,, ^ 


a/3u> — q 


— 1 w —dw w(w — l)(w —-) 


rr)/H = 0, (2.5) 


where /(^p^)$(^pk) = /(s) and q = + q )• Note that the condi¬ 

tion 

^> + |+^i + | + ^ + | = a + /3+l (2.6) 


is satisfied. 
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Equation (l2~5l) with condition (ES|) is called the Heun equation HEU. It has 
four singular points 0, 1, a -1 , oo, all the singular points are regular. The following 
Rienrann’s P-symbol show the exponents. 



Up to here, we have explained how to transform the equation of the Inozemtsev 
model into the Heun equation. Conversely, if a differential equation of second order 
with four regular singular points on a Riemann sphere is given, we can transform 
it into equation 112.51) with condition EH with suitable h ( i = 0,1,2,3) and q 
by changing a variable w —> and a transformation / —> w ai (w — 1 ) a2 (w — 

o-~ 1 )“ 3 /- It is known that if a ^ 0,1 then there exists a solution r to the equation 
a = ei~el ( e * (* = I) 2,3) depend on r). Thus the parameter r is determined. The 
values ei,e 2 ,e 3 and E are determined by turn. Hence we obtain a Hamiltonian of 
BCi Inozemtsev model with an eigenvalue E starting from a differential equation 
of second order with four regular singular points on a Riemann sphere. 


2.2. Trigonometric limit. In this section, we will consider a trigonometric limit 
(r —> y/ —loo). We introduce a parameter p = exp(7T\/—lr), then p —> 0 as r —» 
\/—loo. (Note that the parameter p is different from the one in nao 

If p —> 0, then e\ —> |7r 2 , e 2 —> — )r7r 2 , e 3 —> — and a —> 0, and the relation 
between x and w (see EH ) tends to w = sin 2 ttx as p —> 0. Set 


d 2 t r 2 

Ht = j ^ + Hh +1) - ^— 
dx z sin nx 


+ l\(l\ + 1)- 


t I of /b + f , h + d 

L t — w{w - lH — ^ + (-1 r )~j -b 

t dw z w w — 1 dw 


lo + 




(lo + h + 2) 2 
4w(w — 1) 


}■ 


(2.7) 

( 2 . 8 ) 


Then H —> Ht — + I) an< I equation 12.511 tends to ( Lt — ^ — 

| k(k +1 ))f(w) = 0 as p —> 0. The operator Ht is nothing but the Hamilton¬ 

ian of the BC\ trigonometric Calogero-Moser-Sutherland model, and the equation 
(Lt — C)f(w) = 0 (C is a constant) is a Gauss hypergeometric equation. 

Now we solve a spectral problem for Ht by using hypergeometric functions. We 
divide into four cases, the case l 0 > 0 and 1 1 > 0, the case l 0 > 0 and l\ = 0, 
the case l 0 = 0 and l\ > 0, and the case lo = h = 0. For each case, we set up 
a Hilbert space H, find a dense eigenbasis, and obtain essential selfadjointness of 
the gauge-transformed trigonometric Hamiltonian. The Hilbert space H plays an 
important role to show holomorphy of perturbation in p , which will be discussed in 
section 0 We note that the case lo = 0 and h > 0 comes down to the case Zo > 0 
and h = 0 by setting x —> x + 


2.2.1. The case lo > 0 and l\ > 0. Set 

$(x) = (sin7nc) io+1 (cos7rj:) il+1 , Ht = ^(x)^ 1 Ht$(x), 
then the gauge transformed Hamiltonian Ht is expressed as 


H t = 


= - jl - 2d (i ° + 1) " S ” - (il + 11 ” j ± + (t 0 + i, + 2)V. (2.9) 

ax 1 \ sin ttx cos ttx / ax 
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By a change of variable w = sin 2 7rx, we have 
7~Lt — 7T 2 (2777- -b Iq ~b l\ ~b 2) 2 


= — 47T' 


! |u>( 1 — tw)- 


( 2 . 10 ) 


+ ^-((/o + h + 2) + l)u^ + m(m + lo + h + 2)| 


for each value to. Hence the equation Hr — 7t 2 (2to + Zo + Zi + 2) 2 is transformed 
into a hypergeometric equation. Set 

ipm(x) = ip!^ 0 ’ ll \x) = c m G m (lo + li + 2, 2 ^° 2 + 3 ;sin 2 nx') (to € Z> 0 ), (2.11) 

where the function G m (a, (3',w) =2 -P\(—to, a + to; ft; w) is the Jacobi polynomial 
of degree m and 

I 7t(2to + ?o + l\ + 2)T(to + ^0 + W + 2)r(Zo + to + ^) 

Cm = Y TOir(?n + Zi + |)r(; 0 + |) 2 


is a constant for normalization. Then 


Hripmix) = 7r 2 (2m + lo + h + 2) 2 ip m (x). 


We define the inner products 


( 2 . 12 ) 


(fi9)= dxf(x)g(x), (f,g)<s> = dxf(x)g(x)\$(x)\ 2 . 

Jo Jo 

Then {ip m (x),'ip m '(x))^, = S m ,m'- Set 

H = j/ : R. —> C : measurable with J \f(x)\ 2 \$(x)\ 2 dx < + 00 , 
f(x) = f(x + 1), f(x) = f(-x) a.e. z} 


(2.13) 


(2.14) 


and define an inner product on the Hilbert space H by (•,•)$■ Then the space 
spanned by functions {rp m (x)\m £ Z> 0 } is dense in H. For f(x), g(x) £ HnC 00 (IR), 
we have 

(H T f(x),g(x))<j, = ((H T f{x))®{x),g(x)®{x)) = (H T (f(x)$(x)), g(x)$(x)) 

= ( f(x)$(x),H T (g(x)$(x))) = (f(x),H T g(x))<s>. 

(2.15) 

It follows that the operator Ht is essentially selfadjoint on the space H. 


2.2.2. The case Iq > 0 and li = 0. Set 

$(x) = (sin7ra;) ,0+1 , Hr = $(x)~ 1 Ht$(x), 
then the gauge transformed Hamiltonian is expressed as 




(lo + 1) cos7ra;\ d 


)d^ + {l ° + lf 
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By a change of variable w = sin 2 7 rx, we have a hypergeometric differential equation, 


Ht — it 2 {2m + l Q + l) 2 

= -47t 2 |w(1 - w )-^p + ( ^° 2 + 3 - ((/o + 1) + l)w) ^ +m(m + Z 0 + 1)}- 

(2.16) 

for each m. Set 

V’m(z) = c£C^ +1 (cos7ra;) (m G Z> 0 ), (2.17) 

where the function C^{z) = 2^1 (—m, m + 2v, i/ + ^; 1^) is the Gegenbauer 

polynomial of degree m and = ^/ 2 0+ (™^^2+2)'^° +1 ) 2 ■ Then 

ftTlp%(x) = TT 2 {m + l Q + l)Vm( X )> 

and (^>®(x),^®,(a;))$ = <5 m)m /, where the inner product is defined as 12.ldll for 
>I>(x) = (sin7rx) io+1 . 

There are relations between Gegenbauer polynomials and Jacobi polynomials. 
More precisely, 4>2mi x ) = ipm^^ix) and t/^+i (x) = (cos7r x)4>m'°\x) (m G Z> 0 ). 
Set 

H = |/ : M —> C : measurable with J |/(x)| 2 |3>(x)| 2 dx < + 00 , 

fix) = fix + 2), f{x) = /(-x) a.e. x|, (2.18) 

H+ = {/ G H|/(x) = f{x + 1) a.e. x}, 

H_ = {/ G H|/(x) = -f{x + 1) a.e. x}, 

and inner products on the Hilbert space H and its subspaces H+, H_ are given by 
(•, •)$. Then H + _L H_ and H = H + © H_. 

The space spanned by {t/>®(x) |to G Z> 0 } is dense in H. For /(x), g{x) G H n 
C°°(]R), we have (Wr/(x),j(a;))t = (/(*), Wt 5 (x))$ similarly to 12.1511 . and it 
follows that the operator Ht is essentially selfadjoint on the space H. 

Similar results hold for subspaces H + and H . In fact the space spanned by 
functions {^) m {x)\m G 2Z>o} (resp. {^ m (x)|m G 2Z>o + 1}) is dense in H + (resp. 
H_) and the operator Ht is essentially selfadjoint on the space H + (resp. H_). 


2.2.3. The case Iq = 0 and l\ > 0. Although the case Iq = 0 and l\ > 0 comes down 
to the case Z 0 > 0 and l\ = 0 by setting x —► x + we collect results for the l 0 = 0 
and 1 1 > 0 case for convenience. Set 

$(x) = (cos7rx) il+1 , Ht = ^(x )- 1 Ht ( I>(x), 


then the gauge transformed Hamiltonian is expressed as 


Ht dx 2 2n 


ih + 1) sin 7 txn d ,^22 

)d^ + ih + 1) n ■ 


COS 7 XX 


Now we set 

^mi X ) = Cm C'm +1 (shl7rx) (?7l G Z> 0 ), 
where C^(z) is the Gegenbauer polynomial appeared in (ET7I) and 
~ G ’ _ / 2^i + Km+i 1 + l)m!rp 1 +l )2 , 

c m — Y r(m+2G+2) ’ tIU311 

= n 2 {m + h + 1 fip%{x), 


(2.19) 
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and (' ip ® (x), ^>®,(x))$ = S m ,m', where the inner product is defined as (12.131) . Set 


H = 


{f- 


C : measurable wiht / |/(x)| 2 |<b(x)| 2 dx < +oo, 


fix) = f(x + 2), /( x) = f(-x + 1) a.e. xj, 

H+ = {/£H|/(x) = /(x+l) a.e. x}, 

H_ = {/ £ H| f(x) = -f(x + 1) a.e. x}. 

Here the inner product on the Hilbert space H is given by (•,•)$. Then H + _L HA 
and H = H+ ©H . The space spanned by functions {ip^ (x)\m £ Z>o} is dense 
in H, and the operator TLt is essentially selfadjoint on the space H. 

2.2.4. The case Iq = 0 and l\ = 0. In this case, the trigonometric Hamiltonian is 
Hr = —£s- Set 4>(x) = 1, U T = H t = --£ 7 , ip m {x) = V2cosrmrx, p m (x) = 
■\/2sinTO7rx (to £ Z>i), and ipo(x) = 1. Then 

H T i>m{x) = Tr 2 m 2 ip m {x), (m £ Z> 0 ), 

H T Vm{x) = TT 2 m 2 tp m (x), (to £ Z >0 ), 

('tpm(x),1pm'(x))^ = (ipm{x),<Pm'(x))<s, = 

(lpm(x), Vm'{x))<S> = 0 , 

where the inner product is defined by (tm . Set 


H 


|/ : R —> C : measurable with j |/(x)| 2 dx < + 00 , 
fix) = f{x + 2) a.e. xj, 


Hi = {f(x) £ H 
H 2 = {f(x) £ H 
H 3 = {f{x) £ H 
H 4 = {fix) £ H 


fix) = fix + 1), fix) = fi-x) a.e. x}, 
fix) = fix + 1), fix) = -fi-x) a.e. x}, 
fix) = -fix + 1), fix) = fi-x) a.e. x}, 
fix) = -fix + 1), fix) = -fi-x) a.e. x}. 


Then the spaces H, are pairwise orthogonal and H = ®^ =1 Hj. 

The space spanned by the functions {ip m ix)\m £ Z> 0 } and {p m (x)|?n £ Z> 4 } is 
dense in H and the operator TLt is essentially selfadjoint on the space H, and also 
on subspaces H i; (z = 1, 2, 3,4). 

2.3. Perturbation on parameters a and p{= exp(7r\/— It)). 

As was explained in section I2~121 eigenvalues and eigenfunctions of the Hamilton¬ 
ian H are obtained explicitly for the case p = 0. In this section, we apply a method 
of perturbation and have an algorithm for obtaining eigenvalues and eigenfunctions 
as formal power series in p. 

Since the functions p(x + Wj) (z = 0,1,2, 3) admit expansions H8.16D and pm 
the Hamiltonian H (see (10) 1 admits the expansion 


H = Ht + Ct + ^ Vkix)p k 


( 2 . 20 ) 


fc=1 


where Ht is the Hamiltonian of the trigonometric model defined in (I2.7II . 14(x) 
ik £ Z>i) are even periodic functions with period 1, and Ct = — \ X)i=o + 1) 
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is a constant. Note that for each k the function 14 ( x) is expressed as a finite sum 
of cos 2irnx (n = 0,..., k). 

First we consider the case (fa > 0, l\ > 0), (fa > 0, fa = 0) or (fa = 0, l\ > 0). 
Set 

(1p m (x)$(x) (fa > 0, fa > 0) 

Vm = < 1p%(x)®(x) (fa > 0, fa = 0) 

(x)$(x) (fa = 0,h>0) 

for to £ Z>o- Then v m (to £ Z>o) is a normalized eigenvector of Ht on H. Let 
E m be the eigenvalue of Ht w.r.t. the eigenvector v rn . i.e., 

((2 to + fa + l\ + 2) 2 (fa > 0, l\ > 0) 

E m = \ (tu -Mo + l) 2 (fa > 0, fa = 0) 

((m + l\ + l) 2 (fa = 0, l\ >0). 

Then we have E m ^ E m >, if m ^ m'. 

We will determine eigenvalues E m (p) = E m + CT + J2T=i Em^p k and normalized 
eigenfunctions v m (p) = v m + J^kLi E™' c mjm' v ™'P k °f the operator H T + C T + 
Vk(x)p k as formal power series in p. In other words, we will find E m (p) and 
v m (p) that satisfy equations 

OO 

(H t + C t + ^2 V k (x)p k )v m (p) = Em(p)Vm(p), (2.21) 

fe=l 

(Vm(p),V m (p)) = 1, (2.22) 


as formal power series in p. 

First we calculate coefficients of z2m' rn ,Vm ' = Vk(x)v m • Since Vk{x) is a 
finite sum of cos 2mrx (n = 0,..., k) and the eigenvector v m is essentially a hyper- 

r k\ 

geometric polynomial, coefficients d x m ! rn , are obtained by applying the Pieri formula 
repeatedly. For each to and k, d^ m , ^ 0 for finitely many to'. 

Now we compute Em^ and c^ m , for k > 1. Set c^ m , = By compar¬ 
ing coefficients of v m 'p k , it follows that conditions are equivalent to 

following relations 

/’V' {fc— k'} j{k'} \ _ 1 {k—k’} p{ k'} 

2-^k , — l\2-^m" L 'm,m" %" ,m') Z^k' — l ■ C/m 


„{*} 


Em Em 


k -1 


M 


% = -o EE 


{k '} ffc-fc'} 


k '=1 m" 


k -1 


bW = v V - V MM 

m / ^ / j m,m" m",m / j m,m m 


(m! ^ to) (2.23) 

(2.24) 

(2.25) 


fc'=l 


fe'=i 


Note that the denominator of 12.2311 is non-zero because of non-degeneracy of eigen¬ 
values. Then numbers and Em ^ are determined recursively from (EUl-ESIl 

uniquely. It is shown recursively that for each to and k, #{m" \ c^ m ,, 0} is finite 

and the sums on (E2H1-E2H1) in parameters m" are indeed finite sums. 

Therefore we obtain “eigenvalues” E m (p) and “eigenfunctions” v m (p) of the op¬ 
erator H as formal power series in p. At this stage, convergence is not discussed. 
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Now consider the case lo = li = 0. Though there is degeneracy of eigenvalues 
on the full Hilbert space H, the degeneracy disappears when the action of the 
Hamiltonian is restricted on H, for each i £ {0,1,2,3} and the calculation of 
perturbation works compatibly on each space H,. Hence the calculation is valid for 
the case lo = h = 0 . 


Let us discuss perturbation for equation 12.511 on the parameter 



e 2 ~e 3 A 
ei-e 3 J ' 


We consider the case lo > 0 and l\ > 0. Write 


where 


Lt = I w(w — 1 ) 


3 d 

L = (aw — 1 )L t + aw(w — l)(l 2 + 77 ) - h aqiw + £, 

1 aw 

l i - 9 \ a \ (Jo T l± + 2 ) 2 


dw 2 


lo + i l 


w w — 1 / dw 


Qi = 


(lo + W + I2 ~ I3 + 3 )( 1 q + h + I2 + I3 + 4 ) — (Iq + li + 2) 2 


£ = 


E 


4(ei - e 3 ) 


(Jo + ^2 + 2 ) 2 ^ U.(li + l)e 3 

4 J-—' 


i =0 


4(ei - e 3 ) 


(2.26) 


(2.27) 


Then the equation Lf(w ) = 0 is equivalent to 12.51) . 

We are going to find eigenvalues and eigenfunctions of L as perturbation on a. 
Set 4>m(w) = u m —2 E\ (—m, m + lo + h + 2; 2 i ° 2 + 3 ;u>), where sin 2 nx = w. We will 
use the following relations later which can be found in |5|. 


where 


T 7 / \ ( 2 m + 1 0 + h + 2) 2 ~ ( ^ 

L T iPm(w) = -^- ym(w), 

W1p m (w) = Am^m+^w) + -B m ^ m (w) + 


w(w - 1 )—ijj m (w) = A' m i[ m+ i(w) + B' m ijj m (w) + C' m ip m -i(w). 
dw 


(m + Iq + l\ + 2 )(m + 2 /^_f 3 ) 

(2m + Jo + li + 2)(2m + Jo + li + 3) 
m(m + I 0 + I 1 + 2 )(m + ^j±3) 
(2m + Jo + li + 2)(2m + Jo + li + 3) 
2 m(m + lo T li + 2 ) + 2 } 2 ‘ (Jo + li + 1) 
(2m. + Jo + li T l)(2m + lo + } + 3) 

,/ _ m(m + lo + li + 2 )(Zo — li) 

(2m + lo + li + l)(2m + 10 + 11+3) 
m(m + 2il 2 +1 ) 

( 2 m + lo T li + 2 ) ( 2 m + lo + li + 1 ) 

,, _ m(m + lo + li + 2)(m + 2il 2 +1 ) 

(2m. + Iq + li + 2)(2m + Iq + li + 1) 


<7 _ (2m+/ 0 +/i+2) 2 <7 / \ _ <7 I v^°° 

oet — j , C’mya) — -r a ; 

»m(a) = U m + Efcli £ m , c| n fc} m ,M m 'a fe , and c^jn = 0 for k > 1 and all m. We will 
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determine £m^ and c^ m , (k > 1) to satisfy Lf(w) = 0. Substituting in (12.201) with 
£ = £ m (a), the following relations are shown: 




where 


- ( c * _i_ ~{k — 1} D 

f i'-'to'+I i JDri 


f — f 

i -'m ( ~'m 

£m^ = ^m,m+l^m+1 + <5k~l,oB m + 


+ (m! ^ m), 

(2.28) 
(2.29) 


3 

A-m — (^2 ^m-A-m Ql-^n 


(m + l 0 + h + 2)(m + ^^)(m + 

(2t77 + /q H - H - 2) (2771 + /q H - H - 3) 


— (^2 H - 9 )^ m £"mJ3m Ql 

, 7 , 7 , o\ , (2Z() + 3)(/ 0 + Zl + 1) , 7 , 7 , o\ 

— — y 2 m(m + Iq + + 2 ) H---)(777(777 + to + ti + 2 ) 

(Zo + h + h — h + 3)(Zo + Zi + h + Z 3 + 4)^ 

~r~ ^(2?ti + /o H - Zi ~b 3)( 2777 . + Zo H - Zi -h 1)^ 

(Z 2 H - §)(Zo — Zi)(77l + Iq + l\ -f- 2)777 
(2777 + Zo H - Zi H - 3) (2771 + Zo H - Zi + 1) 

Cm = — (Z 2 + -)C f m — £ m Cm — QlC m 

_ 777(777 + 2 ^hl)( m + Zo+*l-*2+*3 + l )( m + fo+/l-i 2 -/3 ) 

(2m + Zo T li + 2)(2m + + 7i + 1) 

Solving the recursive equations (12.2812.2911 . £ m {a) and u m (a) are obtained as formal 
power series in a. By expanding ei, e?, and e 3 as series in a, eigenvalues of the 
operator L are obtained as formal series in a. 

Now we compare two expansions of eigenvalues (E m (p) <t=> £ m (a)) and eigen¬ 
functions (v m (p) O u m (a))- From the formula a = lGpJX^li ( i+p&»-i ) (see [TSJ 
§21.7]), it follows that a is holomorphic in p near p = 0 and admits an expansion 
a = 16p-f 0{p 2 ). Hence the formal power series u m (a) and £ m (a) are expressed 
as the formal power series in p, and coefficients of p k on u m {a) (resp. £ m (a)) 
are expressed as linear combinations of coefficients of a 1 (l < k) on u m (a) (resp. 
£ m (a)). Set Vm{p) = (1 — aw)~'- l2+1 ^ 2 u m (a) and E m (p) = £ m (a), then they also 
satisfy equation OB as formal power series in p. Since the coefficients are de¬ 
termined by the recursive relations uniquely, it follows that E m (p) = E m (p ) and 
v m (p) = Cm(p)vm(p) , where C m (p) is a formal power series in p. Note that C m (p) 
appears from a difference of the normalization. 

In summary, the perturbation on the variable p is equivalent to the one on the 
variable a. 


3. Perturbation on the L 2 space 
Throughout this section, assume Iq > 0 and l\ > 0. 
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3.1. Holomorphic perturbation. In this subsection, we will use definitions and 
propositions written in Kato’s book J3J freely. The main theorem in this subsection 
is Theorem 13. (il As an application, we show convergence of the formal power series 
of eigenvalues E m (p) (resp. £ m (a)) in p (resp. a) which are calculated by the 
algorithm of perturbation explained in section IQ 

We denote the gauge-transformed Hamiltonian of the BC\ Inozemtsev model by 
H(p), i.e., 


where 


H{p) = <h(x) 1 ° ( — 

{ (sin 7rx) 

(sin 7rx) 

\ 

(COS 7 TX) 

i 


d 2 

dx 2 

3 

+ h(h + 

i=0 

1 )p(x 

+ OJi 

)) 0 

(3.1) 

^0 + 1 j 

(COS 7Tx) Il + 1 

{lo > 

0, /I 

>0); 


^o+l 


( lo > 

0,Zi 

= 0); 


/l +1 


(lo = 

0,Zi 

>0); 




(lo = 

0,Zi 

= 0), 



and p = exp(7T\/~lT). Let I4(x) be the functions in q> , Ct be the constant in 
C2B , and TLt(= d^x) 1 Hj’<f‘(x)) be the gauge-transformed Hamiltonian of the 
BC\ Calogero-Moser-Sutherland model. Then the operator H(p) is expanded as 


H(p) = H T + C T + J2 V k( x )P k - 

i 


Note that H( 0) = Ht + Ct- The functions 14 (x) satisfy the following lemma. 


Lemma 3.1. Let s be a real number satisfying s > 1. Then there exists a constant 
A such that 14 (x) < As k for all k £ Z>i and x £ R. 

Proof. The functions 14 (x) are determined by relations (18.1(11) . Write 

OO 

c(p) = v *p k 

k=l 

oo q 2ti 

= ^ 7r2 (y~!(^o(^o + 1) + h(h + 1))^ P 9r , (^- 2 ) 

n= 1 ^ 

+ (i!(i2 + I) + i3(ia + 1))( n<A tp.)). 

Then the function Y^k= 1 14(2-') is uniformly evaluated by coefficients of majorant 
series C(p) = Y^kLi VkP k for x £ R. Since the convergence radius of (13.211 in p is 1, 
we obtain the lemma. □ 


Proposition 3.2. The operator H(p) (— 1 < p < 1) is essentially selfadjoint on 
the Hilbert space H. 

Proof. Fix the parameter p (—1 < p < 1). Set W(p) = H(p) — H( 0). Then the 
operator W (p) acts as multiplication of a real-holomorphic periodic function of x 
with period 1. We denote the multiplication operator W(p) by W(p\x). 

From El, an inequality W(p\x) < C(p) is shown for all x £ R. Hence the 
boundedness ||IT(p)/|| < C(p)||/|| for f £ H is shown. By definition of the inner 
product on H, the operator W(p) is shown to be symmetric on the whole space H. 
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In sections essential sclfadjointness of the gauge-transformed trigono¬ 

metric Hamiltonian Hr is illustrated. Thus the operator H( 0) is also essentially 
selfadjoint. Therefore essential selfadjointness of the operator H(jp) = H(0) + W(p ) 
is derived from the stability theorem for essential-selfadjoint operators (see 0 V- 
§4.1]). ' ' □ 

Let T denote the closure of a closable operator T. Then H(p) (— 1 < p < 1) 
is the unique extension of H(p) to the selfadjoint operator. Since the symmetric 
operator W(p) is defined on the whole space H, the domain of the operator H(p) 
(— 1 < p < 1) coincides with the one of H( 0). 

Proposition 3.3. The operators H(jp) form a holomorphic family of type (A) for 

-1 < p < 1 . 

Proof. Let s be a real number satisfying s > 1. From Lemma rm we have 
||I4(x)/|| < As fe ||/||, because the measure of an interval (0,1) is 1. On the other 
hand, the operators Vk(x) are defined on the whole space H and symmetric. Hence 
the proposition follows from Kato-Rellich theorem 0 VII-§2.2, Theorem 2.6]. □ 

For the case p = 0, the operator .ff(O) coincides with the gauge-transformed 
Hamiltonian of the trigonometric BC\ Calogero-Moser-Sutherland model up to 
constant. All eigenfunctions of the operator H (0) were obtained explicitly in section 
The spectrum <j(H( 0)) contains only isolated spectra and the multiplicity of 
each spectrum is 1 or 2. Then the resolvent R((,H( 0)) = (( — .ff(O)) -1 is compact 
for C a{H(0)). 

From Theorem 2.4 in 0 VII-§2.1] and Proposition 1301 in this paper, we obtain 
the following statement. 

Proposition 3.4. The operator H(p) has compact resolvent for all p such that 

-1 < p < 1 . 

Let a(H(p)) be the spectrum of the operator H(p). From Theorem 6.29 in 0 
III-§6.8] and Proposition 13.41 in this paper, it follows that 

Proposition 3.5. The spectrum a(H(p)) contains only point spectra and it is dis¬ 
crete. The multiplicity of each eigenvalue is finite. 

Combining Theorem 3.9 in 01 VII-§3.5], Propositions 13.31 and EHI in this paper 
and the selfadjointness of H(p), the following theorem is shown. 

Theorem 3.6. All eigenvalues of the Hamiltonian H(p) (—1 < p < 1) can be rep¬ 
resented as E m (p) ( m £ Z>o), which is real-holomorphic in p £ (—1,1), and E m { 0) 
coincides with the (m + l)sf smallest eigenvalue of the trigonometric Hamiltonian 
H( 0), which was obtained in section 1301 explicitly. The eigenfunction v m (p,x) of 
the eigenvalue E m (p) is holomorphic in p £ (—1,1) as an element in the Hilbert 
space H, and eigenvectors v m (p,x) form a complete orthonormal family on H. 

As an application of the theorem, convergence of the formal power series of 
eigenvalues in the variable p obtained by the algorithm of perturbation is shown. 

Corollary 3.7. Let E m (p) (to £ Z>o) be the eigenvalue of the Hamiltonian H{p) 
of the BC\ Inozemtsev model obtained by the algorithm of perturbation (see section 
IP) . and v m (p) be the eigenvector of the eigenvalue E m {p). If \p\ is sufficiently 
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small then the power series E m (p) converges, and the power series v m (jp) converges 
as an element in the Hilbert space. 

Proof. From Theorem rnn for each m £ Z> 0 and p £ (—1,1), there exists real- 
holomorphic eigenvalues E m (p) and normalized eigenfunctions v m (p,x) £ L 2 which 
converge to the trigonometric ones as p —> 0. Since the eigenvalues E m (p) and the 
eigenfunctions v m (p,x) are holomorphic in p near p = 0, there exists e m £ R.>o 
such that E m (p) and v m (p, x) are expanded as the series in p and they converge on 
|p| < e m . The convergence of v m (p,x) is as an element in H. Eigenvalues E m (p) 
and eigenfunctions v m (p,x) satisfy the following relations for \p\ < e m , 

OO 

( H t + C T + ^2 V k{x)p k ')$(x)v m (p,x) = E m (p)®(x)v m (p,x), 
fc=l 

{®(x)v m (p,x),$(x)v m (p,x)) = 1. 

These equations are the same as remrrm From uniqueness of the coefficients 
obtained by perturbation, it is seen that E m (p) = E m (p) and $>(x)v m (p, x) = v m (p). 
Hence the convergence of E m (p) and v m (p) are shown. □ 

We also obtain the holomorphy of perturbation on the variable a from Theorem 
EH because p is holomorphic in a near a = 0 (see section 

Corollary 3.8. Let£ m (a) (m £ Z>o) (see section Wlft) be the eigenvalue andu m (a ) 
be the eigenvector of the Heun operator L (see IP.PfiV ) obtained by the algorithm for 
perturbation from Lt (see \2.27j ). If |a| is sufficiently small then the power series 
£ m (a) converges, and the power series u m (a) converges as an element in L 2 space. 

3.2. Properties of the eigenvalue. In this subsection, we will show some prop¬ 
erties of eigenvalues. First we discuss the multiplicity of eigenvalues. Under some 
assumptions, it is seen that eigenvalues never stick together. We also introduce 
some inequalities for eigenvalues. 

Theorem 3.9. Assume Iq > 1/2 or l\ > 1/2. Let E m (p) (m £ Z>o) be the 
eigenvalues of H(p) defined in Theorem. \H .61 Then E m (p) E m >(p) (m ^ m!) for 
— 1 < p < 1. In other words, eigenvalues never stick together. 

Proof. Assume E m (po) = E m '(po) for some m < m! and —1 < po < 1. Let us 
consider solutions of a differential equation (H(p 0 ) — E m (p 0 ))g(x) = 0. Exponents 
of solutions at x = 0 (resp. x = 1/2) are Iq + 1 and — Z 0 (resp. Zi + 1 and 
— li). From the condition (Iq > 1/2 or l\ > 1/2), solutions g(x ) satisfying the L 2 
condition f 0 \g(x)\ 2 dx form a one-dimentional space. Therefore the multiplicity of 
the operator H(po) on the Hilbert space H with the eigenvalue E m (po) is one. 

From Proposition 13.51 the spectrum a(H(p 0 )) is discrete. Hence there exists 
ei £ R>o such that H(po) has exactly one eigenvalue in the interval (E m (po) — 
ei,E m (po) + ei). Since the operators H(p) (—1 < p < 1) form a holomorphic 
family, there exists £2 £ R>o such that H(p) has exactly one eigenvalue in the 
interval (E m (jPa) - ei/2, E m (po ) + ei/2) for p such that p 0 - e 2 < p < po + £2 (see 
(U V-§4.3, VII-§3.1]). Set f(p) = E m (p) — E m >(p). Then the function f(p) is real- 
holomorphic in p £ (—1,1), f(po) = 0, and /( 0) 0. From the identity theorem for 
the real-holomorphic function f(p), it follows that the point p = po is an isolated 
zero, i.e., there exists e £ R>o such that E m (p) E m ’(p) f° r 0 < \p — pq\ < e. 
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Hence if |p — po|(^ 0) is sufficiently small then values E m (p), E m i(p ) belong to 
the interval ( E m (p 0 ) — ei/2, E m (p 0 ) + e\/2) and E m {p) ^ E m i ( p ). This shows that 
H(p ) has no less than two eigenvalues in the interval {E m (po) — ei,E m (po) + ei), 
and it contradicts. Therefore, we obtain the theorem. □ 

Corollary 3.10. Assume Iq > 1/2 or l\ > 1/2, then E m (p) < E m '(p) for m < m! 
and —1 < p < 1. 

Proof. From the labelling of the eigenvalues E m (p ), it follows that E m { 0) < E m i (0) 
for to < m!. It is seen in Theorem Em that values E m (p) and E m '(p) never stick 
together for — 1 < p < 1. Hence E m (p) < E m '(p) for m < m' and — 1 < p < 1. □ 

Now we introduce another inequality. 

Theorem 3.11. Assume l 2 = h = 0. Let E m (p) (to £ Z>o) he the eigenvalues 
of H(p) defined in Theorem, Vi.(A For each to, the eigenvalue E m (p) is monotonely 
increasing for 0 < p < 1 and monotonely decreasing for — 1 < p < 0. 

Proof. Set H(p) = H( 0) + YlkLi ^2 k{x)p 2k ■ From the formulas 18.161 l8J71) . it is 
seen that 14 k{x) > 0 for all k £ Z >0 . Hence ^ (EfcLi ^2 k{x)p 2k ) > 0 for 0 < p < 1 
and | (E^i V 2k (x)p 2k ) < 0 for -1 < p < 0. 

Let v m (p,x) be the normalized eigenfunction of H{p) labelled in Theorem 13.61 
By definition we have H(p)v m (p,x ) = E m (p)v m (p, x) and 

(H(p)v m (p,x),v m (p,x)}^ = E m (p). (3.3) 

By the way, we have (f^v rn (p,x),v m (p,x))$, = 0, which is obtained by differentiat¬ 
ing the equality (u m (p, x), v m (p, x))<p = 1 in p. 

Let us differentiate the equality (TOll . From the right hand side, we have 
-f^E m (p). From the left hand side, we obtain 

d 00 d 

( (^ XJ ^ 2fc (x)p 2k ) Vm ( P, x) , V m (p, *) ) $ + 2 (H(p)v m (p, x),—V m (p, x)) $ 

/>1 7 OO J 

= J ^YlV 2 k{x)p 2k ^\v m {p,x)\ 2 \^(x)\ 2 dx + 2E m (p)(v m (p,x), — v m (p, x))^ 

rl 7 00 

= J £ VMP 2k ) \*rn(p, x)\ 2 Mx)\ 2 dx. 

Here we used the selfadjointness of H(jp). Thus 

7 /»1 j OO 

^£^m(p) = J ^^V 2 k(x)p 2k ^\Vrn{p,x)\ 2 \<f>(x)\ 2 dx. 

If 0 < p < 1 then the integrand is nonnegative for all x £ [0,1]. Hence ^ E m {p) > 0. 
The inequality E m {p ) < 0 for — 1 < p < 0 is shown similarly. □ 

4. Perturbation on the space of holomorphic functions 

In the previous section, we obtain eigenfunctions v m (p, x) of the BCi dimensional 
Inozemtsev model as elements in the Hilbert space. On the other hand, for the case 
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Zoj 6 Z>o and h = h = 0 , holomorophy of the function v m {p,x) in the variable 
x for sufficintly small |p| is shown by applying the Bethe Ansatz method Ca¬ 
in this section, we show holomorphy of the function v m ( p , x) in the variable x 
for sufficintly small \p\ for the case Z 0 , l\ G R>o, h, Z 3 £ R. The proof of holomorphy 
in x is similar to the one for the elliptic Calogero-Moser-Sutherland model of type 
An, which was done in [5j. The paper [Hj would be helpful to understand contents 
in this section. 

Let v m be a normalized eigenfunction of the trigonometric Hamiltonian 44(0) 
and E m be the corresponding eigenvalue which were obtained explicitly in sections 
12.2.112.2.41 

The action of the resolvent (44(0) — £) _1 on the Hilbert space H is written as 
(44(0) £) 'y ] c m v m = y ] (44 m £) c m v m , 

m m 

for J2m CmVm e H. Hence if C ^ er(44(0)), then the operator (44(0) — () _1 is 
bounded. Set W(p) = 44 (p) — 44(0). Then the operator W(p) is bounded and 
||lF(p)|| —► 0 as p —> 0. If ||(44(0) — C) _ 1 4L(p)|| < 1 then the resolvent (44(p) — ^) _1 
is expanded as 

{H{p) - C )” 1 = (1 + (44(0) - C)- 1 VF(p))” 1 (44(0) - C )” 1 

00 • 14 1) 

= £(-(£(°) - C)- 1 W(p)) J (44(0) - C)- 1 - 

4=0 

Let T be a circle which does not bump the set cr(44(0)) and let r = dist(T, er(44(0))). 
Then there exists po > 0 such that ||W(p)|| < r and T does not bump the set 
cr{H{p)) for all \p\ <po- Set 

Then the operator Pr{p) is a projection to the space of linear combinations of 
eigenfunctions in the Hilbert space H whose eigenvalues are inside the circle T. Fix 
E m € cr(44(0)). Since the set er(44(0)) is discrete, we can choose a circle T m and 
Pm G R>o such that T m contains only one element E m of the set cr(44(0)) inside 
T m and the projection P m (p) = Pr m {p ) satisfies \\P m (p) -P m ( 0 )|| < 1 for \p\ < p m . 

For the case (Zo > 0, h > 0), (Zo > 0, Zi = 0) or (Zo = 0, Zi > 0), multiplicity 
of every eigenvalue of the operator 44(0) on the Hilbert space H is one. Then a 
function P m {p)v m for each m is an eigenvector of the operator 44 (p) and admits an 
expression 

P m {p)v m = C(p,m)V m (p , x), (4.2) 

( c (p,m) is a constant) for sufficiently small |jp|, because the operator 44 (p) preserves 
the space P m {p )H and v m (p, x) is an eigenvector of 44 (p) with the eigenvalue E m {p) 
(see Theorem 13.011 . 

For the case lo = h = 0, the Hilbert space H is decomposed into Hi ® H 2 ® 
H 3 ® H 4 (see section 12. 2. 411 . For each space, there is no degeneracy of eigenvalues 
for the operator 44(0). Hence the expression (14.211 is also valid. 

As for the expansion of P m (p)v m in v m i {m! E Z>o), the following proposition 
which is analogous to |5; Proposition 5.11] is shown. 






16 


KOUICHI TAKEMURA 


Proposition 4.1. Let \p\ < 1. Write Pm{jp)v m = Y m ' ■ Bor eac h 171 

and C G K>i, there exists C' G K>o and p * such that coefficients s m ,m' satisfy 
< C"(C|p|) |m_m,|/2 for all m! G Z> 0 and p (|p| < p*). 

We prove this proposition in section [S] To obtain holomorphy of the function 
Pm(p)vm ? we need the following proposition, which we will prove in section |5J 

Proposition 4.2. Let , if m (x)(= u m ) &e i/ie (m + l)st normalized eigenfunction of 
the trigonometric gauge-transformed Hamiltonian 'Hr ■ Let f(x) = Ym=o c miftm(x) 
be a function satisfying \c m \ < AR m (Vm G Z>o) for some A G R>o and R G (0,1). 
If r' satisfies 0 < r' < -Mog-^, then the power series ^m=o c m’4’m{A) converges 
uniformly absolutely inside a zone —P < Qx < r', where Qx is an imaginary part, 
of the complex number x. 

Remark. By the relation w = sin 2 nx, a zone |3x| < log A i n the variable x 
corresponds to a domain inside an ellipse in the variable w with fori at w = 0,1 
and a major axis r = R+R —. 

Holomorphy of the power series v m (p) in variables x and p is shown by applying 
Propositions 14.11 and 14.21 

Theorem 4.3. Let v m (p) be the eigenfunction of the Hamiltonian H{p) obtained 
by the algorithm of perturbation (see section HOI) . Then the power series v m (p) 
is an analytic function in variables x and p for \p\: sufficiently small. More 

precisely for each r G R>o there exists po G R>o such that the series v m (p) converges 
absolutely uniformly for (p,x) G [—po,Po] x B r , where B r = {x G C|| 3 x| < r}. 

Proof. We prove for the case l 0 > 0 and f > 0. Write P m {p)v m = Y m ' s m,m' v m ' 
and apply Proposition 14.11 for the case (7=2. Then there exists C' G M.>o and 
p* such that coefficients s m>m / satisfy < C"(2|p|)l m-m I/ 2 for all p, m' such 

that |p | < p*. 

Let po = min(p*/2, exp(-57rr)/2). Then we have r < ^ log < b log 

and there exists A G R>o such that |s mj771 '| < A{\J 2jp[) m for |p| < po(< p*) 
and m! G Z>q. By Proposition 14.21 the series Y m ' s m,m'V m ' converges absolutely 
uniformly for (p, x) G [—po,Po] xB r . Since the function v m (p ) constructed in section 
12.31 coincides with a function P m (p)u m up to constant and the function P m (p)u m = 
Ym' s m,m'Vm' converges absolutely uniformly for (p, x) G [—Po,Po] x B r: we obtain 
absolutely uniformly convergence of v rn (p). 

For the case (Iq > 0 and l\ = 0) or (l 0 = 0 and l± > 0) or (l 0 = l\ = 0), we 
can prepare alternative propositions to Propositions 14 . ll and 14. 2l and can prove the 
theorem similarly. □ 


5. Algebraic eigenfunctions 

5.1. Invariant subspaces of doubly periodic functions. If the coupling con¬ 
stants Zoi ^i) ^ 2 , ^3 satisfy some equation, the Hamiltonian H (see ED) preserves a 
finite dimensional space of doubly periodic functions. In this section, we look into a 
condition for existence of the finite dimensional invariant space of doubly periodic 
functions with respect to the action of the Hamiltonian H (see Proposition 15. 111 . 
After that, we investigate relationship between invariant spaces of doubly periodic 
functions and L 2 spaces. 
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Proposition 5.1. Let o.i be a number such that di = —U or cti = L + l for each i £ 
{0,1,2,3}. Assume - X)Lo “*/ 2 e z >o set d = - Xu=o c^/ 2 - ^4o,ai,a 2 ,a 3 

&e Z/ie d + 1 dimensional space spanned by j ( a <f(x) ) ai ( ) a2 () a3 p(z) n I n = 

0,..., dj , where a(x) is the Weierstrass sigma-function defined in and at (x) 

(i = 1,2, 3) are the co-sigma functions defined in Q. Then the Hamiltonian H 
(see \2.H) ) preserves the space Va 0 ,ai,& 2 ,a 3 - 


Proof. Set z = p(x), <&(.z) = (z —e\) ai ^ 2 (z —e 2 ) a2 ^ 2 {z — ef) a3 ^ 2 , and H = $(z) 1 o 
H o <j>(z). Then 


TT A ! \ / \ l , f ^ V-' Oi{ T 2 d 

H = -A(z- ei )(z - e 2 )(z - e 3 )< — + — “ 


dz ■ 


+ { - (di + d 2 + d 3 — l 0 ) (di + d 2 + d 3 + lo + 1) z 
+ ei(d 2 + d 3 ) 2 + e 2 (di + d 3 ) 2 + e 3 (di + d 2 ) 2 }. 


Let Vd +1 be the space of polynomials in z with degree at most d. From formulas 
633 , it is enough to show that the operator H preserves the space Vd+i. 

The action of the Hamiltonian is written as 

H(z - e 2 ) r = -4^(r + 7i)(r + 7 2 )(2 - e 2 ) r+1 

+ ((e 2 — e 3 )(r + d 2 + di)r + (e 2 - ei)(r + d 2 + d 3 )r + q)(z — e 2 ) r 

+ r(r + a 2 - i)(e 2 - e 3 )(e 2 - ei )(z - e 2 ) r_1 ), 

(5.1) 

where q' = (ei(d 2 + d 3 ) 2 + e 2 (di + d 3 ) 2 + e 3 (di + d 2 ) 2 ) + 627172 , 

71 = (di + 012 + d 3 — Zo)/2 and 72 = (di + a .2 + d 3 + Iq + l)/2. Hence the operator 
H preserves the space of polynomials in z. Since do = —lo or Zo + 1, it follows 
that (r + 7i)(r + 72) = 0 for r = d, (d = — X^=o cfei/2). Hence the coefficient of 
(z — e 2 ) r+1 on the right hand side of EH> is zero for the case r = d. Therefore the 
operator H preserves the space Vd+ 1 - □ 


Proposition 5.2. With the notation in Proposition ] 5. l\ assume d = — Oj/2 £ 
Z>o- For the trigonometric case ( p = 0), the eigenvalues of the Hamiltonian H 
on the finite dimensional space Va Q ,ai, 5 2) a 3 are written as {n 2 (2r + do + di) 2 — 

— Si=0 4" l)}r=0,...,d- 

Proof. As p —> 0, we have e\ —► 27t 2 /3, e 2 —► —7 t 2 /3, and e 3 —> —7 t 2 /3. In this 
case the coefficient of (z — e 2 ) r ~ 1 on the right hand side of 115.111 is zero for all r. 
Hence the operator H acts triangularly. Then the eigenvalues appear on diagonal 
elements. By a straightforward calculation, the coefficient of (z — e 2 ) r on the right 
hand side of <E3J is written as 

2 3 

7r 2 (2r + a .2 + d 3 ) 2 —— ^ ' li(h + 1). 

i=0 

Here we used relations d 2 — on = k(k + 1) (z = 0,1,2,3). Therefore the eigenvalues 
of the Hamiltonian H on the finite dimensional space V r 5 0 , 5 1 ,a 2 , 5 3 are {7r 2 (2r + 
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d 2 + d 3 ) 2 - ^-Yj\=oh{h + l)}r=o,...,d- By replacing r —■> d - r, we obtain the 
proposition. □ 

5.2. Algebraic eigenfunctions on the Hilbert space. In this subsection we 
investigate a relationship between the Hilbert space H and the finite dimensional 
space Vs 0 , 53 ,s 2 , 5 3 which was defined in the previous subsection. Throughout this 
subsection, assume Iq > 0 and fa > 0 . 

Let us consider the case di £ {—fa, k + 1} (i = 0,1, 2, 3) and d = — Xt=o € 
Z>o- It is easily confirmed that if do > 0 and di > 0 then any function f(x) in 
Va 0 ,a lt & 2 ,&3 is square-integrable on the interval ( 0 , 1 ). 

Set T4 0 ,5i,a 2 ,a 3 = (4^y| f( x ) e v a o,ai,fi 2 ,s 3 }> where $(x) is the ground state of 
the trigonometric model which was defined in section l2~2l 

Proposition 5.3. Let oti £ {—fa, fa + 1} for i = 0,1, 2, 3. If do > 0, di > 0 and 
d = — E?=odj/2 £ Z>o, then V& 0 , 53 , 53,53 C H. For the case (Iq > 0 and fa = 0J 
or (Iq = 0 and 1 1 > 0) we have Va 0 , 53 , 52,53 C for i = + or — , and for the case 
Iq = h = 0 we have Va 0l 53 ,a 2 , 5 3 C for i = 1, 2 , 3 or 4. 

Proof. We consider four cases, i.e., the case Iq > 0, h > 0, the case Iq > 0, l\ = 0 
case, the case Iq = 0 , l\> 0 , and the case Iq = l\= 0 . 

First, we prove the proposition for the case Iq > 0, l\ > 0. The numbers do and 
di must be chosen as do = Zo + 1 an d di = /i + l from the condition do, di, Zo, h > 0 . 
Now we check that, if f(x) £ V 5 0 , 53 , 5 2 ,s 3 , then f{x) satisfies the definition of the 
Hilbert space H (see 12.14H - ). Square-integrability of the function f(x) follows from 
the condition do > 0 and di > 0. Periodicity and symmetry of f(x) follow from 
the condition — Xu-o cfei/2 £ Z. Hence f(x) £ H. 

For the case Iq > 0 and l\ = 0, do and di are chosen as do = Zo + 1 an d (di = 0 
or 1). Then V 5 0 , 5 i, 5 2 , 5 3 C for i = + or — is shown similarly. Note that the sign 
of i is determined by whether functions in V 5 0 , 53 , 5 2 , 5 3 are periodic or antiperiodic 
(for details see section 15. 2.211 . 

For the other cases, the proofs are similar. □ 

Remark. A function in V 5 0 , 53 , 5 2 , 5 3 is multi-valued in general, and so we should 
specify branches of the function. In our case, the analytic continuation of the 
function near the real line should be performed along paths passing through the 
upper half plane. 

In the finite dimensional space V 5 0 , 5 i, 5 2 ,a 3 ) eigenvalues are calculated by solv¬ 
ing a characteristic equation, which is an algebraic equation, and eigenfunctions 
are obtained by solving linear equations. In this sense, eigenvalues in the finite 
dimensional space are “algebraic”. 

Now we figure out properties of the spaces V 5 0 , 53 , 53 , 5,3 contained in the Hilbert 
space H. We divide into four cases. 

5.2.1. The case Iq > 0 and l\ > 0. Let d 2 £ {—fa, fa + 1} and a 3 £ {—fa,fa + 1}. 
From Proposition 15.31 if d = —(Iq + l\)/2 — 1 — (d 2 + a 3 )/2 £ Z> 0 then the d + 1 
dimensional vector space Vz 0 +i,q+i, 5 2 , 5 3 is a subspace of the Hilbert space H. We 
will show that the set of eigenvalues of the gauge-transformed Hamiltonian H(p) 
(see (13.Ill ) on the space V/ 0+ i,/ 3 +i,< 5 2 ,a 3 is the set of small eigenvalues of H(p) on 
the Hilbert space H from the bottom. 
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Lemma 5.4. For the trigonometric case p = exp(7r\/"^Tr) =0, if d = —(Iq + 
l\)/2 — 1 — &2 — <53 £ Z >0 ZZien Z/ie seZ of eigenvalues of the gauge-transformed 
Hamiltonian H( 0) on the finite dimensional space Vz 0 -|_i,Zi+i,a 2 ,a 3 | p= o coincides with 
the set of small eigenvalues of H{ 0) on the Hilbert space H from the bottom. In 
other words, the m-th smallest eigenvalue of H(ff) on H is also an eigenvalue on 
Vi o+ i th+lt&2t&3 \ p = 0 , if and only if m < d + 1. 

Proof. From Proposition 15.21 eigenvalues of the trigonometric Hamiltonian H on 
the finite dimensional space V&o,&i,& 2 ,c «3 are written as {tt 2 ( 2 r + Iq + h + 2) 2 — 
Xu =0 h(h + l)} r= o,...,d- On the other hand, from equality (12.121) and the limit 
H —> Ht — + 1) as p —> 0, eigenvalues of the gauge-transformed 

trigonometric Hamiltonian H(ff) = | p=0 on the Hilbert space H are 

written as {7r 2 (2 r + l 0 + h + 2) 2 — J2i=o ^*(^» + l)}re z> 0 - Therefore the lemma 
follows. □ 

From the previous lemma, eigenvalues of the trigonometric gauge-transformed 
Hamiltonian H{ 0) on the finite dimensional space are {i? r (0) | r = 

0 ,... ,d}, where the values E r (p) are defined in Theorem Id. (il It is obvious that the 
eigenvalues of the operator H(p) on the finite dimensional space Vz 0 +i i i 1 +i i 5 2i a 3 are 
continuous in p. Hence if — 1 < p < 1 then the eigenvalues of the gauge-transformed 
Hamiltonian H(p) on the finite dimensional space V r ; 0 +i,; 1 +i,a 2 ,a 3 coincide with 
{E r (p) | r = 0,..., d}. By applying Corollary Id.101 for the case Zq > 1/2 or the case 
Zi > 1/2, we obtain the following theorem. 

Theorem 5.5. Assume d = — (Iq + h)/2 — 1 — (&2 + &z)/2 € Z>o- If flo > 1/2 
and l\ > 0/ or (Iq > 0 and l\ > 1/2/, then the set of eigenvalues of the gauge- 
transformed Hamiltonian H(p) on the finite dimensional space Hz 0 +i ! ; 1 _|_i i , 5 2i a 3 co¬ 
incides with the set of small eigenvalues of H{p) on the Hilbert space H from the 
bottom. In other words, the m-th smallest eigenvalue of H(p) on H is also an 
eigenvalue on Vi 0 +i y i 1 +i,a 2 ,a 3 , if and only if m < d+ 1. 

5.2.2. The case Iq > 0 and l\ = 0. Let € {—I 2 , h + 1} and 013 £ {— Z3, 13 + 1}. If 
d = —(Zo + l)/2— (a 2 + <T 3)/2 £ Z>o, then the d+ 1 dimensional space Vz 0 +i i o, 5 2 ,a 3 
is a subspace of the space H + , and if d = — (l 0 + 2)/2 — (5 2 + ds)/2 £ Z> 0 , then 
the d + 1 dimensional space Vz 0 +i,i,a 2 ,a 3 is a subspace of the space H_. 

Similarly to Theorem ESI the following statement are shown. 

• If d = —(Zo + l)/2 — ( 0:2 + a 3)/2 £ Z>o, Iq > 1/2, and — 1 < p < 1, then the set 
of eigenvalues of the gauge-transformed Hamiltonian H(p) (see (Id.Ill ) on the finite 
dimensional space Vz 0 +i,o,a 2 ,a 3 is the set of small eigenvalues of H(p) on the space 
H + from the bottom. 

• If d = — (Zo + 2)/2— (d 2 + Q!3)/2 £ Z>o, Iq > 1/2, and — 1 < p < 1, then the set of 
eigenvalues of H{p) on the finite dimensional space I/ 0 +i,i,a 2 ,a 3 is the set of small 
eigenvalues of H(p) on the Hilbert space H_ from the bottom. 

5.2.3. The case Zo = 0 and l\ > 0. Let a .2 £ {—Z 2 , h + 1} and 013 £ {—Z 3 ,13 + 1}. If 
d = -(Zi+l)/ 2 -(a 2 -l-a 3)/2 £ Z> 0 , then Vo,ii-(-i, 5 2 ,a 3 C H + and dim Vb.Zi+i.as,^ = 
d + 1. If d = -(Zi + 2)/2 - (d 2 + & 3)/2 £ Z> 0 , then Vi,z 1 +i,a 2 ,a 3 C H_ and 
dim V\ z 1+ i a 2 &3 = d + 1. We can confirm similar statements to the ones in section 

Em' 
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5.2.4. The case Iq = 0 and h = 0. Let a 2 £ {—h,h + 1} and (S 3 £ {—^ 3,^3 + 1}. 
If d = —{a .2 + ck 3)/2 £ Z>o, then Vb,o,a 2 ,a 3 C H 1? C H 2 , dim Vb,o,a 2 ,a 3 = 

d+1, and dim Vi ) i i a 2i 5 3 = d. If d = —1/2—(a 2 +a: 3)/2 £ Z>o, then Vi i o,5 2) a 3 C H 4 , 


Vo,i,a 2 ,a 3 C H 3 , and dim Vi,o,a 2 ,a 3 = dim Vo,i,a 2 ,a 3 = d + 1- 

5.3. The case of nonnegative integral coupling constants. If the coupling 
constants lo,h,l 2 , h are nonnegative integers, the model satisfies some special prop¬ 
erties. Specifically, it admits the Bethe Ansatz method m and the potential has 
the finite-gap property mmm- In this subsection, we reproduce several results 
more explicitly for the case lo,li,h,h are nonnegative integers. Note that some 
results were obtained in m. 

Throughout this subsection, assume l t £ Z> 0 for * = 0,1,2,3. Assume (3 i £ Z 
(i = 0,1, 2, 3) and — JA_ 0 A/2 £ Z> 0 . Let V^ o ^ ^ ^ be a vector space spanned 


^ {(^) /3l (w) /32 (^) /33 ^) rl }„=o,...,-ELoft/2 End We S6t VfoAM = 


e where $(z) was defined in section^ Let a* £ {-l it k + 


1} (i = 0,1,2,3) and 


U & 


Oq,Oi ,02,03 


U & 


oq,o:i,0:2,013 


^OQ ,di ,02 ,03 5 


Vx- 

,{ 0 }, 


Oq , 1 — Ol ,1 — 02,1 — O3 5 


Vqcq ,01,02 ,03 5 


^1- 

,{ 0 }, 


Oq , 1 — Oi ,1 — 02,1 — O3 5 


Xa=0 “*/2 G ^<0; 
Z)Lo “i / 2 e z> 2 ; 

otherwise, 

Xa=o “j / 2 e ^<0; 
Xa=o “*/ 2 e 
otherwise. 


If Iq + l\ + I 2 + h is even, then the Hamiltonian H (see (12.Ill 1 preserves the spaces 
bb- lo,—li,—I 2 ,—I 31 bb- Iq,— £,i 2 +i 4 3 +ij 1^—L4 3 +i( bb_ io4i+i,i 2 -l-i,—Z 3 • The gauge- 
transformed Hamiltonian H(p) (see (13.Ill ) preserves the spaces C/_j 0i _j lt _j 2 _j 3 , 

L r -i 0 ,-z ll i 2 +i,z 3 +i, I/-io,h+i,-L43+i, ^-L,h+i,i 2 +i,-i 3 - H b) + h + h + I 3 is odd, 
then the Hamiltonian H preserves the spaces |_i, [/_; 0 i _z 1) z 2 _|_i,-z 3 , 

C/_; 0j i 1 -)_i j _; 2i _; 3 , [/; 0 +i ! __; 1 ,_i 2i _; 3 , and the gauge-transformed Hamiltonian f?(p) 
preserves the spaces C/_; 0 _ ili _; 2i z 3+ i, U-i 0 - ht i 2+ i-i 3 , U-i 0th+1 -i 2 -i 3 , and 
b/o + L —A — ^2.-^3 ■ 

Now we present results on the inclusion of Ua 0 ,ai,a 2 ,a 3 in a Hilbert space. We 
consider eight cases. 


5.3.1. The case Iq > 0, l\ > 0 and Iq + h + l 2 + I 3 is even. If l 2 + I 3 — lo ~ h > 2, 
then V; 0+ i ) ; 1+ i i _; 2) _i 3 C H and dim V; 0+ i i i 1+ i i _i 2i _ i3 = ( l 2 + I 3 - lo ~ h)/2. 

5.3.2. The case Iq > 0, l\ > 0 and Iq + h + l 2 + I3 is odd. If l 2 — Iq — h — I3 > 3, 
then L/Q+ijZi+i,—; 2 ,; 3 +i C H and dim Vz 0 + i,z 1 +i,-z 2 ,z 3 +i = ( l 2 — Iq — h — I3 — l)/2. 
If h-h-h-h > 3, then Vz 0 + 1 ,z 1+1 ,; 2+1 ,-z 3 C H and dim V io+1M+1 ^ 2+1 _* 3 = 
(I3 ~ lo ~ h ~ h ~ l)/ 2 - 

5.3.3. The case Iq > 0 , l\ = 0 and I0 + I1 + I2 + I3 is even. If I2 + I3 — Iq > 2 , then the 
space V; 0 +i i i i _j 2 j _i 3 is a subspace of the space H_(c H) and dim VJ 0 -|_i,i,-z 2 ,-z 3 = 
{I2+I3—I0)/ 2 . H 1 2 —Io—Iq > 2 , then Vj 0 +i i o,-z 2 ) z 3 +i C H + and dim Vj 0 + i j o,-z 2 ) z 3 +i = 
(h—lo—h)/^- IH3—I0—I2 > 2 , then Vj 0 +i i o,z 2 +i,-z 3 C H+ and dim Vz 0 +i,o,z 2 +i,-z 3 = 
(I3 — lo ~ k)/2. 
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5 . 3 . 4 . The case Iq > 0 , l\ = 0 and Iq + h + h + I3 is odd. If I2 + h — lo > 1 , then 

Vio+i,o,-i 2 -h C H + and dim Vz 0 +i,o,-z 2 ,-z 3 = ( I2 + I3 — I0 + l)/ 2 . If I2 — I0 — I3 > 3 , 
then hz 0 +i,i,-z 2) ; 3 +i C H_ and dim Vz 0 +i,i,-z 2 ,z 3 +i = ( I2 — — 1 )/ 2 . If I3 — I0 — 

I2 > 3 , then Vz 0 +i,i,z 2 +i,-z 3 C H_ and dim V) 0 +i,i,z 2 +i,-z 3 = (I3 — Iq — I2 — l)/ 2 . 

5 . 3 . 5 . The case Iq = 0 , h > 0 and Iq + h + I2 + I3 is even. If I2 + I3 — h > 2, then 
Vi,Zi+i,-z 2 ,-z 3 C H_ and dim V\,i 1 + i-i 2 -i 3 = (Z2 + Z3-Zi)/ 2 . UI2 — I1 — I3 > 2 , then 
bb,Zi+i,-Z 2 ,i 3 +i C H + and dim Vo,Zi+i,-z 2 ,z 3 +i = {l2 — h — h)/ 2 - If h ~ h — h > 2 , 
then Vo.Zi+i,z 2 +i,-z 3 C H + and dim Vo,Zi+i,z 2 +i,-z 3 = {I3 — h ~ h)/ 2 . 

5.3.6. The case Iq = 0, Zi > 0 and Iq + h + I 2 + h is odd. If I 2 + I 3 — h > 1, then 

bb,Zi+i,-z 2 ,-z 3 C H + and dim Vo^+i-i^-^ = (Z 2 +13 ~ h + l)/2. If I 2 — li — I 3 > 3, 
then ,-; 2 ,; 3 -|-i C H_ and dim V’i i ; 1 -|_i i _; 2) ; 3 _|_i = (I 2 — Zj — Z 3 — l)/2. If I 3 — l\ — 

I 2 > 3, then t^i,; 1 +i,; 2 +i ,-; 3 C H and dim Vi^+i^+i.-Zs = (Z 3 — Zi — Z 2 — l)/2. 

5.3.7. The case lo = h = 0 and Iq + h + h + I 3 is even. In this case, the spaces 

Vo,o-i 2 ,-i 3 , b),i,—z 2 ,—z 3 ? 

{ Vi,o,-z 2> ; S+l) Vo,i ,-W 3 + 1 (h > I 3 ) 

hi,o,z 2 +i,-z 3I Vo,i,/ 2 +i,-z 3 {I 2 < h) 

{0} (h = k) 

are subspaces of the Hilbert space H. 

5.3.8. The case Iq = l± = 0 and Iq + h + h + I 3 is odd. In this case, the spaces 


^1,0, — 12, —13 5 

12,—13 > 




1 V\p-i 2 j 3+ \, 

I/o.o ,-z 2 ,z 3 +i 

(h > I 3 ) 


Vi,l,j 2 + l,-l 3 , 

Vb,o,z 2 +i,-z 3 

{h < I 3 ) 


l {0} 


(l 2 = h) 


are subspaces of the Hilbert space H. 


6. Examples 

In this section, we show examples that illustrate results of this paper. 

6.1. The case Iq = 1, li = 2, Z 2 = 0, Z3 = 8. In this case, the Hamiltonian is 

d 2 

H = —+ 2p(ai) + 6p(x + 1/2) + 72 p(x + t/2). 

Set $(3:) = (sin7nc) 2 (cos7nr) 3 and consider the gauge-transformed operator H(p) = 
< h(a;)^ 1 o H o $(2:), where p = exp(ny/—lr). By the trigonometric limit p —> 0, we 
obtain H( 0) = TLt — ^f-, where 

d 2 ,2 cos 7 tx 3sin7nr N d 9 , ,, 

= -)—+ 257T 2 . 6.1 

dx* sin ttx cos ttx dx 

The Hilbert space H and the inner product are defined by (12.13112.1411 . Then 
the eigenfunction of the gauge-transformed trigonometric Hamiltonian H( 0) with 
the eigenvalue (2m + 5) 2 7r 2 — is the Jacobi polynomial = ipm’ 2 \x) (see 

(12.1111 ). The functions (ipm{x))m£Z >0 form an orthogonal system on the Hilbert 
space H. Let H(p) be the selfadjoint extension of H(p). From Theorem 13.61 and 
Corollarv l3.10l it is shown that all eigenvalues of H (p) (—1 < p < 1) on the Hilbert 
space H can be represented as E m (p) (m £ Z>o), which is real-holomorphic in 
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p £ (— 1 , 1 ) and E m ( 0 ) = ( 2 m + 5) 2 7r 2 — Moreover, E m (p ) < E m '(p) for 

m < ml and — 1 < p < 1. 

Define spaces of doubly periodic functions as follows: 



and set : f(x) £ ^Then the Hamiltonian H 

preserves the spaces 14 , 3 , 1 ,-8, H- 1 ,- 2 , 1 ,- 8 , 14,- 2 , 0 ,-8, F- 1 , 3 , 0,-8 and the gauge- 
transformed Hamiltonian H(p) preserves the spaces 14, 3 , 1 ,- 8 ; V- 1 ,- 2 , 1 ,- 8 , 14,- 2 , 0 ,- 8 , 
h— 1 , 3 , 0 ,—8- Among the spaces 14, 3 , 1 ,- 8 ; fo- 1 ,- 2 , 1 ,- 8 , 14,- 2 , 0 ,- 8 ; 1 , 3 , 0 ,—8> only 

the space 14, 3 , 1,-8 is a subspace of the Hilbert space H. The eigenvalues of 
the gauge-transformed Hamiltonian H(p) on the space 14 , 3 , 1,-8 are written as 
llei — 9e2 ± 2-\/l06ef + 73eie2 + 46e|, where e\ = p(l/2) and e 2 = p(( 1 + r)/2). 
From Theorem 15.51 we obtain Eq{p) = llei — 9e2 — 2y / 106e 2 + 73eie2 + 46e| and 
E-[ ip) = llei — 9e2 + 2-\/l06e 2 + 73eie2 + 46e|. Hence the smallest eigenvalue and 
the second smallest one on the Hilbert space H are obtained algebraically. 


6.2. The case Iq = 1, Zi = 2, 4 = 1, I 3 = 0. In this case, the Hamiltonian is 

d 2 

H = — 2 -f- 2p(x) -f 6 p(x + 1/2) + 2 p(x + (t -f-1)/2). 

Set $(x) = (sin7rx) 2 (cos7rx) 3 and consider the gauge-transformation H(p) = <f>(x) _1 o 
H o <f>(x). The Hilbert space H is defined similarly to the case Iq = 1. 1-\ = 2 , Z 2 = 

0 ,13 = 8. Define 



and set ^ f(x) £ Then the Hamiltonian H pre¬ 

serves the spaces y_i,_ 2 ,-i,0) V-i,-2,2,11 14,-2,-1,1 and the gauge-transformed 
Hamiltonian H(p) preserves the spaces y_i,_2,-i,0) H-1,-2,2,1, 14,-2 ,-i,i- But none 
of spaces V_i,_2,-i,o, lAi-2,2,1, 14 ,-2,-1,1 is included in the Hilbert space H. 

6 . 3 . The case Iq = 1 , 4 = 0 ,12 = 4 , I3 = 1 . In this case, the Hamiltonian is 

rp 

H = — ——— + 2p(x) + 20p(x + (t + l)/2) + 2p(x -1- t/2). 
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Set $( 2 ;) = sin7ra and consider the gauge-transformation H(p) = <I>(x) _1 o Ho^(x), 
where p = exp(7iV— It). By the trigonometric limit p —> 0, we obtain H( 0) = 
Ht — 87t 2 , where 


Ht 


dx 2 


— 27r 


2 cos 7 tx 3 sin 7rx 


sm nx 


cos nx 


d 2 
— + 47T. 
dx 


( 6 . 2 ) 


The Hilbert space H, its subspaces H + , H_, and the inner product are defined 
by (12.131 12.1811 . Then the eigenfunction of the gauge-transformed trigonometric 
Hamiltonian H( 0) with the eigenvalue (m + 2) 2 7r 2 — 87r 2 is written by use of the 
Gegenbauer polynomial ipm( x ) ( see 12.17t h The functions (i/>^(:r)) m6 z >0 (resp. 
(ipm( x ))m£ 2 Z> 0 , (ipm( x ))m£ 2 Z> 0 +i) form an orthogonal system on the Hilbert space 
H (resp. H+, H_). Let H(p) be the selfadjoint extension of H(p). From Theorem 
EZD and Corollary mm it is shown that all eigenvalues of H(p) (—1 < p < 1) on 
the Hilbert space H can be represented as E m (p), which is real-holomorphic in 
p € (—1,1) and E m (0) = (to + 2 ) 2 7 t 2 — 87r 2 . Define 


w / ,cr 2 (x)^-4,a 3 (x)^-i 

v -x,o,-4^x = {( TU rr) (-^r) P{x) 


n=0,l,2,3 


v 2 . 


/^cr 1 (x) w a 2 (x)^-4 r a 3 (x)^-i 
^ p(x) , n =o,i 




,_ 4 , 2 =«)’>, 
\ v a{x) (J{x) 

, CTl (x) S , (J 2 (x) , -4 / cr 3 (x) X 2 


T T I \ -4 (WA\*) N .2 . , f 


71=0,1 


and set Vs 


= : f(x) £ v p 0 ,p u p 2 ,p 3 }- Then the Hamiltonian H pre- 


_ f( x ) 

serves spaces V r _i i o,-4,-i, V2 ,i,-4,-i, ^,0,-4,2, R-1.1,-4,2, the gauge-transformed 
Hamiltonian H{p) preserves spaces V-i t o,-4,-i, V2.i,-4,-i, V2,o,-4, 2, H-1,1,-4,2, and 
we have V^o,—4,2 C H+, and V 2 p : — 4,—1 C HT. 

The eigenvalue of the gauge-transformed Hamiltonian H(p) on the space V 2 ,o,- 4,2 
is —12ei — 12 e 2 and the eigenvalues of H(p) on the space V^, 1 ,_ 4 ,_ 1 are written as 
12ei — 3e 2 ± 2y / 39ef^P3eie 2 _ ~6ef. From the results in section Hi. 2. 21 the smallest 
eigenvalue on H + is — 12ei — 12e2 and the smallest two eigenvalues on H_ are 12ei — 
3e 2 ± 2^/39e 2 + 3eie 2 — 6e|. Hence E 0 (p) = —12ei — 12e2, Ei (p) = 12ei — 3e 2 — 
2y / 39ef"+"3e)l( 2 ^ — i6ef, and E 3 (p) = 12ei — 3e 2 + 2-\/ / 39ef^i^3eie 2 — ^6e|. In other 
words, the smallest eigenvalue, the second smallest one, and the fourth smallest one 
on the Hilbert space H are obtained algebraically. 


7. Concluding remarks 

In the paper pQ, Gomez-Ullate, Gonzalez-Lopez, and Rodriguez found square- 
integrable finite dimensional invariant spaces for the BCn Inozemtsev model with 
some coupling constants, and they observed numerically for an example of the BC 2 
Inozemtsev model that the set of eigenvalues on a finite dimensional invariant space 
would coincide with the set of small eigenvalues of the Hamiltonian on L 2 space 
from the bottom. In the present paper, we justified this phenomena concretely 
for the BC\ Inozemtsev model in Theorem 1531 and sections EH 15331 and we 
illustrated plainly with examples in section[5] It would be possible to obtain similar 
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results for the BCn Inozemtsev model by combining the Kato-Rcllich theory and 
the trigonometric limit. 

We comment on relationship to the known facts for the Heun equation. In 
books mini, notions “Heun polynomial” and “Heun function” are introduced. 
Roughly speaking, the Heun polynomial corresponds to a doubly periodic eigen¬ 
function of the BC\ Inozemtsev model, and the Heun function corresponds to a 
square-integrable eigenfunction or a doubly periodic eigenfunction of the BC\ In¬ 
ozemtsev model. We hope that both the Heun equation and the Inozemtsev model 
are elusidated substantially in near future. 


8. Appendix 


8.1. Proof of Pronosit,ion 14.11 In this subsection we prove Proposition 

Lemma 8.1. Let ifmix) be the normalized eigenbasis of the gauge-transformed 
trigonometric Hamiltonian Ht roi defined in 1071) . and write Vk{x)p k = 
H{p)~H{ 0). Then 

4>m(x)V k (x) = E Cm’ll>m'(x), (8.1) 

\m'—m\<k 

for some constants c m ' (| m! — m\ < k). 


Proof. From the Pieri-type formula for the Jacobi polynomials, it follows that 

il>m(x) cos 2nx = c_V>m-i(:E) + c 0 ip m (x) + c + i/j m+1 (x), (8.2) 

for some constants c_, Cq and c+. On the other hand, the function V k ( x ) is a 
polynomial in the variable cos27ra; with degree k. By applying llO) repeatedly, we 
obtain ED- ' □ 

Proposition 8.2. Let \p\ < 1 and QT)j£L 1 Vk(x)p k )ip m (x) — Em'e Z> 0 tm.m'lfm'(x). 
For each C such that C > 1 and C\p\ < 1, there exists C" £ M>o such that 
\im,m'\ <C ,, (C\p\) i2SL =r i±1 . 


Proof. Since the normalized Jacobi polynomials form a complete orthonormal sys¬ 
tem, it follows that f m>m / = (lp m >(x), QlfcLl Vk(x)p k )lfm(x))^. 

If k < | m! — m |, then (Vk{x)p k ipm(x), 4>m'{x)}$ = 0 by Lemma l%Jl and orthogo¬ 
nality. Therefore, 


<Vm'(z), 

V /c=i 1 


OO 

Vw(z),( E Vk(x)p k ^l/Jm{x) 


k=\m' —m\ 


l oo 


E Vk(x)p k 1 pm(x)lp m '(x)$(x)‘ 2 dx 

k=\m' —m\ 


< sup 
xe[o,i] 


OO „1 

y: V k (x)p k / | '1/Jm(x)^m'(x)®{x) 2 \dl 

, i Jo 


k=\m' — m\ 


E Vk \p\ k - 

k'>\m' — m\ 


< 
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where 14 is defined by 13.211 . For the case m! = to, we obtain \t m ,m\ < Efc>i Vfc |p| fc - 
Since the convergence radius of the series Yh n ^nP n is equal to 1, for each C 
such that C > 1 and C\p\ < 1 there exists C" £ R>o such that \tm,m'\ < 
C"{C\p\)\ m '- m \ (to 7 ^ to') and < C ,,, (C'|p|) < C"(C\p\) x / 2 . From an in¬ 
equality (C|p|)l m,_m l < (Cjp|) J for | m' — m\ £ Z>i, we obtain the proposi¬ 
tion. □ 


Proposition 8.3. Let D be a positive number. Suppose dist(£, a(H ( 0 ))) > D. 
Writ e ( H(p ) - C,)~ l ifm{x ) = Y^ m ' t m,m'il’m'{x), where ( H(p ) - C) _1 defined in 
14 - 1 ]) . For each m £ Z>o and C £ R>i, f/iere exists C' £ M>o trod po £ R>o which 
do not depend on £ (but depend on D) such that t mtm i satisfy 

|i m y|<C'(C|p|) |m '- ra|/2 , (8.3) 

for all p, m! such that \p\ < po and m' £ Z>q. 


Proof. Recall that the operator (H(p)—Q 1 is defined hv Neumann series 14.11) . Fix 
a number D(£ R>o) and set X := (£ — -ff(0)) _1 (EfcL;L Vk{x)p k )- From expansion 
14. Ill , there exists pi £ M>o such that an inequality ||X|| < 1/2 holds for all p and 
f such that |p| < pi and dist(£, > D. Then we have (Ei^o -X'*)(^ r (0) — 

C)” 1 = (H(p) - C)" 1 - K we write Em' c m'Vw(z) = (#(0) - C) _1 E m ' c m ' 4 >m’(x), 
then 14,1 < £> _1 |cm'| for each to'. Write Xip m (x) = E m ' im,m'ifm'{x). By com¬ 
bining with Proposition IQ we obtain that for each C such that C > 1 and 
E/Ipi < 1, there exists C" £ R>o which does not depend on £ (but depend on D) 
such that < C"(^±I|p|)(l m '- m l +1 )/2 for \ p \ < pi . 

To obtain Proposition IQ we use the method of majorants. For this purpose, 
we introduce symbols e m (to £ Z) to avoid inaccuracies and apply the method of 
majorants for formal series Em'6Z Cm ' em ' • Bor formal series, we define the partial 
ordering < by the following rule: 


V c (1) e ,< V c (2) e , 

/ , L m'^rn’ 

m'EZ m' EZ 


o Vro', |c£)| < |c^|. 


„( 2 )| 


We will later consider the case that each coefficient (z = 1,2, to' £ Z) is ex¬ 
pressed as an infinite sum. If absolute convergence of is shown for each m! , 
then absolute convergence of e^j for each m! is shown by the majorant. Set Xe m = 
Em' where coefficients 4,m' are defined by Xip m (x) = E m / im.m'Vw (%)■ 

Our goal is to show 18.31) for t m ,m' such that E m ' tm,m'tpm'(x) = E^o 

Since (.ff(0)-C) _: V m (a:) = (-E m -C) _1 4W and |(£ , m -C) _1 | < -D -1 , 
it is enough to show that there exist C* £ R>o and po which do not depend on ( (but 
depend on D) such that t* m m , are well-defined by EiEo -4 e m = E m '£z4,m' e m' 
and satisfy 

(C - 1-1 \ l m '- m l/ 2 

14,m'I <C*(—^|p|) , (8.4) 

for all p and to' such that |p| < po and to' £ Z>q. Set 


Z em, — 




•2m, m'^m' — 


E 

m'EZ 


c" 


4 + 1 y 

\— p ) 


Then we have Xe m <Ze m . Let k £ Z>i. If coefficients of Z k e m with respect to 
the basis {e m ,} converge absolutely, then X k e m is well-defined and X k e m <Z k e m . 
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From the equality Z k e m = EkCO,..., ! ,o-i) z m,A 1 )Z v m, v w ■ ■ • and the 


property z m>TO / = it follows that 

OO OO 1 /* 

E z " em = E E j (E z ^ s ") 

1--1 blm'ci v -Msl-i „ cZ 


gm _ m '_! dg 


= E E 


^ez 


1 


27lV 1 -/|s| = 


f ( V C" , p (H+1) s l/ ') k s m ~ m '~ 1 dse r 

^l =1 1 


k= 1 ra'GZ 

= ^ ^ Zm,m' ? 
ra'GZ 

where p = (^-5-^p) 1 / 2 and 

2ttv / ^1 £ N i (1 -ps)(l -ps” 1 ) - C"'(p — p 3 ) ’ 


2 

Zm,m' 


C"(p — p 3 )s 


;3\ c m—m! — 1 


ds 


(8.5) 


Note that we used a formula Enez9^ +la; " = (l-ga-K 1 --?^" 1 ) ' Equality EH make 
sense for p < p 2 , where P2 is a positive number satisfying inequalities p 2 < 1, 
C (i_p7)^ 21 < 1 and C"p2 < 1. Therefore each coefficient of Ea*Li ^ fce m with respect 
ot the basis {e m >} converges absolutely. It follows that 


^ ' 27 G m ^G m -p ^ ' Z G m f^G m + ^ ( Z m rn fG m '. (8.6) 

k =0 fc=l m' EZ 

Let s(p) be the solution of an equation (1—ps)(l —ps” 1 ) — C"(p —p 3 ) = 0 on the 
variable s satisfying |s(p)| < 1. Then s(p) is holomorphic in p near 0 and admits 
an expansion s(p) = p + c 2 p 2 + ..., and 

1 r C"(i,-f)s^4s (8.7) 


(2tt\/^I) J| s |=i (1 -ps)(l ~ps 3 ) - C"(p — P 3 ) 

where /(p) is a holomorphic function defined near p = 0. Note that relation (18.71) 
is shown by calculating the residue around s = s(p) for the case n > 0 and we need 
to change a variable s —> s _1 and calculate the residue around s = s(p) for the case 
n < 0. By combining 18.5118.71 there exists a positive number such that 


y X k e m <e m + y p/(p)s(p) |m m| e m /. 


( 8 . 8 ) 


/c=0 


m/EZ 


for |p| < p 3 . By combining (18.81) . a relation p = (E±ip) 2 , an inequality E±I < C 
and an expansion s(p) = p + c 2 p 2 + ..., we obtain 18.41) and the proposition. □ 


The following proposition is essentially the same as Proposition 14.11 

Proposition 8.4. Let E m £ cr(H( 0)) and T m be a circle which contains only one 
element E m of the seta(H( 0)) inside it. Letil> m (x) be the corresponding normalized 
eigenfunction. Set P m (p) = — 27T ^_ 1 f r ( H(p ) — ()~ 1 df and write Pmip^mix) = 
Em' s m,m l ' l Pm'[x). For each C £ R>i, there exist C' £ M>o and p* £ M>o such 
that Sm,m' satisfy 

ISm.m'l < C"(C|p|)l m '- m l/ 2 , 
for all p and m' such that |p| < p* and m! £ Z> 0 . 


(8.9) 
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Proof. Since the spectrum cr(H(0)) is discrete, there exists a positive number 
D such that dist(£, a(H(0))) > D for ( £ T m . Write {H(p) — C)~ 1 ip m (x) = 
E m '^ m i m '(0' t l’m , (x). From Proposition 18.81 it is shown that for each C £ R>i 
there exists C* £ R>o and p* £ R>o which do not depend on £(e r m ) such that 

tm,m' satisfy |i m>m '(C)| < C'*(C'|p|) J 2 for all p, m! such that \p\ < p* and 
to' £ Z> 0 . 

Let L be a length of the circle T m and write — 2n ]_ 1 f r (H(p) — ()~ 1 d(if m (x) = 
Era' s m,m'4>m'(x). By integrating the function E m ' *m,m'(C)Vw(z) over the circle 

r m in the variable £, it follows that |s m>m /| < ^(7*(C|p|) 2 for all p, to' such 

that \p\ < p* and rrt! £ Z>g. Therefore, Proposition l8.4l is proved. □ 


8.2. Proof of Proposition 14.21 We prove Proposition 14.21 

Proposition 8.5 (Proposition Let if m (x) be the (m + l)st normalized eigen¬ 
function of the trigonometric gauge-transformed Hamiltonian Hr- Let f(x) = 
E°^L=o c m^m{x) be a function satisfying |c m | < AR m (Vto £ Z>o) for some A £ 
K>o and R £ (0,1). If r' satisfies 0 < r' < log A , then the power series 
0 C m'4 , m{x) converges uniformly absolutely inside a zone —r' < 'Ax < r', where 
Ax is an imaginary part of the complex number x. 


Proof. We prove for the case Iq > 0 and l\ > 0. For the other cases, they are proved 
similarly. We introduce a Rodrigues-type formula for the Jacobi polynomials 

p m { W ) = 1 _ ( w lo+m+1 /2 (1 _ w) h+m +iM _ 

Then Pm (sin 2 7 nc) = where d m = ^ From 

the Stirling’s formula, we have (d m ) 1/,m —> 1 as to —> 00 . Hence it is sufficient 
to show that the power series Em=o c mPm (sin 2 7 ra;) converges uniformly absolutely 
inside the zone —r' < Ax < r'. 

The generating function of the Jacobi polynomials p m (w ) is written as 


Y = 


o{ l+5+S\ ; o+l/2 j l-£+S\h+l/2 ’ 

m =0 M” 2 ~ ) (— ) 

where S = i/(l + £) 2 - 4 £ w. Now we set Em=o Qm{y)^ m = 

Em=o9m(y)C" = and 

OO 

Y Prniy)^ = 


( 8 . 10 ) 


x/(i-€) 2 -4y€ ’ 


1 


m =0 


(VU - 0 2 - (1 - e + V(! - 0 2 - 4 y^yo+h+i ' 


( 8 . 11 ) 


Then it is shown that, if a > 0, then q m {y ) and {y) are polynomials in y of degree 
m with nonnegative coefficients. Hence p m (y) is also a polynomial in y of degree m 
with nonnegative coefficients. Set p m {y) = EfcL 0 Pm^y fc and p m {y) = EEo P'm y k ■ 
From formulas (IHUni EHH> and the nonnegativity, we obtain |p„ | < Prn for all TO 
and k. From the inequality | sin 2 7 rx| < | (-—-)“| for — r' < Ax < r ', it is seen 


2 
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that 


|pm(sin 2 7rx)| < Y (sm 2 TTx) k \ < Y?m I sin 2 nx\ k 
k =0 k—0 


< 




( 8 . 12 ) 


for all m G Z>o and x such that —r' < Ax < r'. 

On the other hand, the series Y^=oPrn({^-—^——) 2 )£, rn 1 with respect to the 
variable £, has radius of convergence e~ 27Tr ; because the singular point of the 
right hand side of iETTTIi which is closest to the origin is located on the circle 
|£| = e~ 2vr . Let r" be a positive number such that r' < r" < ^ log(l /R). 

Then there exists A! G R. such that p m ((— —-§ -) 2 ) < A' e 2 ' nr m . Hence we 

have \ Cm \\p m ({^+^) 2 )\ < AA r (Re 2nr ") m . Since Re 2 * r " < 1, the series 

Ylm=o\ c m\\pm ({-— + 2 e * ) 2 ) | converges. From inequality (18.121) , uniformly abso¬ 
lute convergence of Y^m=o c mPm( sin 2 ttx) inside the zone —r' < ’Ax < r' is obtained. 
Therefore, the proof is complete. □ 


8.3. We note definitions and formulas of elliptic functions. Let ui\ and u > 3 be 
complex numbers such that the value w 3 /aq is an element of the upper half plane. 
The Weierstrass p-function, the Weierstrass sigma-function and the Weierstrass 
zeta-function are defined as follows: 
p(x) = p(x\2ui,2u> 3 ) 

= L + Y ( _I_I_ ) (8.13) 

12 { m,n)^\mO)} y(X ~ 2mUJl ~ 2nUJ3)2 (2mWl + 2na;3 ) 2; ’ 


a(x) =x ( 

(m,n)GZxZ\{(0,0)} 

X 


1 - 


X 


x exp 


( 2 muq 


2 towi + 2nu>3 J 
x 2 

-(- 2 nu >3 ' 2(2muji -(- 2 ?ro ; 3 ) 2 
a'(x) 


C( x ) = 


a{x) ' 


Setting u> 2 = —toi — u> 3 , = p(u>i) and rji = C(w;) for i = 1,2,3 yields the relations 

ei + e 2 + e 3 = r?i + 772 + r? 3 = 0, 

p{x) = ~C(x), (p'(x)) 2 = 4 (p(x) - ei)(p(x) - e 2 )(p(x) - e 3 ), 

p(x + 2cui) = p(x), C(x + 2u>i) = £(x) + 2rji (i = 1,2,3), 
p"(x) 1 t 1 1 


(p'(x)) 2 


A 


x — e\ x — e 2 x — e% 


k)- 


/ . \ _ . ki &i')ki Ci") ,. _ 0 

p(x + Vi) — e, -f- — - (* — 1,2,3), 

p{x) - a 

where i ', i" G {1,2,3} with i' < i", i ^ i' and i 7 ^ i". 

The co-sigma functions <Ji(x) (i = 1,2, 3) are 

a i(x) = exp(—r]ix)a(x + w;)/cr(wi). 


( 8 . 14 ) 
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and satisfy 

(~Tt) 2 = P( x ) ~ e i> (* = 1,2,3). (8.15) 

a(x) 

Set u>i = 1/2, W 3 = r/2 and p = exp( 7 r\/—lr). The expansion of the Weierstrass 
p-function in the variable p is written as 


p(x) 



3 


n —1 


np , . 

--— (cos2n7nr — 1 . 

i - p 2n v ' 


(8.16) 


By setting x —> x+1/2, x + r/2, x+(1 + t)/ 2, the following expansions are obtained 


/ jl \ /i 71 o \ Tip 

p ( x+ 2^ = cosher) _ y - 87r Y 

v 7 n=l 


-((— 1 )" cos2nitx — 1 ), 


p(* + J) = - g7r2 Y 


np 


p(a 


1 + b 


n—1 

oo 




,, cos 27r?ra; — p n 

1 -p 2n ’ 

(—l) n cos 2nnx — p n 


(8.17) 


n=1 


i -p 


,2n 
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